Iiep-lat/9305010| RIKEN-AF-NP-154 



Phases of the three-state Potts model in three spatial dimensions 



Shigemi Ohta 

The Institute of Physical and Chemical Research (RIKEN) 
Wako-shi, Saitama 351-01, Japan 



Abstract 

The three-state Potts model is numerically investigated on three-dimensional simple cubic 
lattices of up to 128'^ volume, concentrating on the neighborhood of the first-order phase tran- 
sition separating the ordered and disordered phases. In both phases clusters of like spins are 
observed with irregular boundaries. In the ordered phase the two different non-favored spins 
are found to attract each other with a long but finite range. As a result, the neighborhoods of 
the non-favored spins are interpreted as domains of the disordered phase. This explains why 
the first-order phase transitions associated with the global symmetry, including the SU(3) 
pure-gauge one, are so weak. 



Finite-temperature systems with global 
Z3 symmetry in three spatial dimensions 
have in common a weak first-order phase 
transition separating ordered and disordered 
phases. Lattice numerical calculations of 
the SU{3) pure-gauge quantum chromody- 
namics (QCD) at finite temperature IH re- 
vealed that a first-order phase transition 
separates a low-temperature color-confining 
phase and a high-temperature non-confining 
phase. The confining phase is disordered 
with the Polyakov line order parameter 
taking zero expectation value. The non- 
confining phase is ordered and the Polyakov 
line takes a finite expectation value falling on 
one of the three Z3 axes in the complex plane. 
The transition is considered weak because the 
latent heat is small and the Polyakov line cor- 
relation length grows. Numerical simulations 
of the three-state Potts model showed 
there is a first-order phase transition separat- 
ing a low-temperature ordered phase and a 
high-temperature disordered one. This tran- 
sition is weak in the sense the latent heat is 
small. However, it is not yet understood why 
these phase transitions are so weak. This let- 
ter reports the author's investigation of this 
problem. 

We consider the simplest form of the three- 
state Potts model on three-dimensional sim- 
ple cubic lattices. At each site, x, of the 
lattice, a spin, ax, is defined. It takes one 
of the three possible states in the Z3 group, 
Co: G {0,1,2}. The thermodynamics of the 
system is described by the partition function, 

Z = exp(-iJ/r) = exp(J 6{ax,ay)), 

where the sum runs the nearest-neighbor 
pairs {xy) and Kronecker's symbol is defined 
as 6{ax,cry) = 1 \i a.^ = (Jy and other- 
wise. The model is invariant under global Z3 
transformations. Here we consider only the 
"ferromagnetic" couplings, J > 0. The effec- 
tive interactions among the Polyakov lines in 
the finite-temperature pure-gauge QCD are 
of this type Q . 

Earlier numerical simulations showed that 
a weak first-order phase transition that sep- 
arates a low-temperature ordered phase and 
a high-temperature disordered phase at J ~ 
0.5505. The phase transition is considered 
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Figure 1: Distribution of the dimensionless 
internal energy density from heatbath simu- 
lations at J = 0.5505 on 32^ (32) and 128^ 
(128) volumes. 

weak because the dimensionless internal en- 
ergy density, e = {5{ax,cry)), shows only a 
small gap at the phase transition, from about 
0.58 to 0.53, compared with the maximum 
possible gap from 1 to 1/3. See Fig. || for 
distribution of this quantity. These results 
were, however, obtained on relatively small 
volumes such as 32^. The two broad peaks 
marked as "32" in Fig. |l] shows the distri- 
bution from such a small- volume simulation 
reproduced by the author. Although the two- 
peak structure is clear enough for establish- 
ing the first-order nature of this phase transi- 
tion, the significant overlapping of the broad 
tails in the middle refiects the fact that flip- 
fiop transitions between the two phases occur 
too frequently and it is hard to tell in which of 
the two phases the system resides at a given 
instant. Because of this difficulty, the natures 
of the individual phases could not be studied, 
nor the reason why the phase transition is so 
weak. 

To avoid this difficulty one has to simu- 
late the system on larger volumes. The au- 
thor wrote a heatbath simulation code for the 
model that runs on an experimental parallel 
computer, APIOOO, built by Fujitsu Labora- 
tory Q . On the largest existing configuration 
of the computer with 16 x 32 = 512 micro- 
computer "cells," the code performs one mil- 
lion heatbath updates of a 128^ volume in 
about 32 hours. Typically a couple of million 
heatbath updates are made at four values of 
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the coupling, J = 0.55025, 0.5505, 0.55075, 
and 0.551. At each of these coupUng values, 
two independent simulations were made; one 
starting from a completely ordered configura- 
tion and the other starting from a completely 
disordered one. 

At the highest temperature, J = 0.55025, 
the ordered-start simulation quickly con- 
verged to the disordered phase within lOOK 
heatbath updates, while the disordered-start 
simulation remained in the disordered phase 
indefinitely through the duration of the sim- 
ulation for more than one million updates. 
In contrast, at the lowest temperature, J = 

0. 551, the disordered-start simulation quickly 
converged to the ordered phase while the 
ordered-start one remained in the ordered 
phase indefinitely. Coexistence of the two 
phases is observed at both of the intermediate 
temperatures, J = 0.5505 and 0.55075: The 
ordered-start simulations remained in the or- 
dered phase and the disordered-start ones re- 
mained in the disordered phase through their 
durations of more than two million heatbath 
updates. The corresponding distribution of 
the internal energy density is shown in Fig. 

1, marked as "128." The two peaks are now 
completely separated, and there is no doubt 
about in which of the phases the system re- 
sides: We are now able to investigate the na- 
tures of the individual phases. 

The two-point correlation function, Cjj(r), 
is defined as the probability to find a pair 
of spin i and j separated by the distance r. 
They are calculated once in every lOOK heat- 
bath updates, and are proven free of auto- 
correlation at least in the relevant range of 
r < 8. 

In the disordered phase the correlations are 
expected to decay exponentially in Yukawa 
form 



exp( 



as the distance r increases toward infinity 
1^]. The constant term piPj is given by 
the product of probabilities, pi and pj, to 
find a site with spin i and j respectively. 
In the disordered phase they should all ap- 
proach 1/3 as the system size increases, be- 
cause the ^3 symmetry is preserved. For 
the same reason, there are only two differ- 
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Figure 2: Correlations in the disordered 
phase at J = 0.5505. The constant terms 
PiPj are subtracted. The curves are obtained 
by the least-x'^ fittings to the Yukawa form 
summarized in Table |l|. 



ent sets of correlations: the diagonal ones, 
Coo = Cii = C22, and the off-diagonal ones, 
Cqi = C02 = C12. Moreover, because of the 
conservation of probability, 

Cioir) + Ca{r) + Ci2{r) = 1, (1) 

for any i, all the correlation mass must agree 
with each other and the amplitudes obey a 
simple relation qoq = an = 022 = — 2aoi = 
— 2ao2 = —2ai2- We numerically confirmed 
that all the correlations are fitted well by 
the Yukawa form, as can be seen in Fig. § 
where the correlations are shown with the 
constant terms subtracted. Within statisti- 
cal errors, the fitting parameters (see Table 
|l|) satisfy the above requirements. The posi- 
tive diagonal amplitudes, an > 0, shows the 
like spins attract each other, in accordance 
with the ferromagnetic interaction. The cor- 
responding correlation mass of ma ~ 0.15 
suggests there exist clusters of like spins of 



Table 1: Two-point correlations in the disor- 
dered phase at J = 0.5505. Two-parameter 
least-x^ fit to the Yukawa form for the range 
r > 3. 
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Figure 3: Diagonal correlations in the or- 
dered phase at J = 0.5505. means the fa- 
vored spin and 1 and 2 are non-favored. The 
constant terms are subtracted. The curves 
are from least-x'^ fitting summarized in Ta- 
ble |. 

the size of several lattice spacings. Indeed 
in the spin distributions one sees many such 
clusters. It is interesting to note that these 
clusters are neither smooth nor convex in 
shape, but are complex and concave. 

In the ordered phase there is one favored 
spin which we label as 0, and two non-favored 
ones which we label as 1 and 2. No ana- 
lytic result is known about the behavior of 
the correlations. The current numerical re- 
sults show the Yukawa form fits them well 
also in this case. See Fig. ^ and 0, where the 
constant terms are subtracted again. The fit- 
ting parameters are given in Table ^. The 
positive amplitudes of the diagonal correla- 
tions, an > 0, are in accordance with the fer- 
romagnetic interaction. Smaller correlation 
mass for the favored spin, 0, compared with 



Table 2: Two-point correlations in the or- 
dered phase at J = 0.5505. Two-parameter 
least-x^ fit to the Yukawa form for the range 
r > 3. The favored spin is 0. 
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Figure 4: Off-diagonal correlations in the or- 
dered phase at J = 0.5505. The constant 
terms are subtracted. The curves are from 
least-x'^ fitting summarized in Table ^. 

those of non- favored spins, 1 and 2, may not 
be statistically significant, but nonetheless is 
consistent with the fact that the favored spin 
dominates the volume in the ordered phase. 
At this couphng, J = 0.5505, about 60 % 
of the volume is covered by the favored spin, 
while the remaining is evenly split between 
the two non-favored ones. On the other hand 
the negative amplitudes of the correlations 
between the favored spin and either of the 
two non-favored spins, aoi < and ao2 < 0, 
means the favored spin repel the non-favored 
spins. These suggest that there form clus- 
ters of like spins, which are indeed observed 
in the spin distributions. Again the clusters 
have irregular boundaries |^. 

The correlation between the two different 
non-favored spins, Ci2{r), is the most inter- 
esting. It starts from zero at the origin as it 
should be, but then overshoots the asymp- 
totic value piP2, and approaches it from 
above (see Fig. ^.) This means an attrac- 
tive force acts between the two different non- 
favored spins. The attractive part alone can 
be fitted by the Yukawa form if we limit the 
range of fitting to r > 3 (Fig. §.) Note that 
the correlation mass, mi2 = 0.03 =b 0.03, is 
much smaller than those for the other cases, 
the smallest of which is moo = 0.14 ± 0.02. 
It is not appropriate, however, to consider 
this as a truly long-range correlation with 
zero correlation mass. Again because of the 
conservation of probability (||) the correlation 
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mass mi2 in the long-distance limit r — > cx) 
must be equal to the smaller of the masses 
moi and rrin. Much smaller value of mi2 we 
are seeing must come from subtle balance of 
the two terms Coi(r) and Cii(r) in the range 
we are looking at, and has to give way to the 
true correlation mass at some longer range. 
For this reason, we will call this attractive 
part as the "middle-range" attraction in the 
remainder of this letter. 

Yet it is useful to explore what would hap- 
pen if we had a really long-range attraction 
between the two different non- favored spins: 
such an attraction is likely to cause insta- 
bility of the ordered phase, because the re- 
pulsion against the non-favored spins from 
the favored ones can then be compensated 
by paring the two different non-favored spins 
attracting each other. Hence we expect that 
the smaller the correlation mass mi2, the 
weaker the phase transition. Indeed in the 
spin distribution of this phase one finds many 
clusters of the non-favored spins. And such 
a non-favored cluster almost always accom- 
pany clusters of the other non-favored spin 
in its neighborhood. It should be noted that 
such pairing of the clusters of the two differ- 
ent non-favored spins is a way to maintain the 
global Z3 symmetry which requires the vol- 
umes occupied by the two non-favored spins 
must be equal. 

As was mentioned earlier, these clusters 
of the non- favored spins appear in irregular 
shapes. The neighborhood of such a clus- 
ter of complex-shaped clusters of the non- 
favored spins can be interpreted as an island 
of the disordered phase in the sea of the or- 
dered phase. This is clearly seen by looking 
at the local-averaged internal energy density, 
ei{r), which is defined in the sphere of radius 
r with spin i at its center. See Fig. || for its 
dependence on the radius. In the disordered 
phase, the quantity does not show any de- 
pendence and stays at its asymptotic value of 
~ 0.53. In the ordered phase, the one around 
the favored spin approaches the asymptotic 
value of ~ 0.58 monotonously from above, 
while the other around the two non-favored 
spins starts from a value smaller than in the 
disordered phase and approaches the asymp- 
totic value monotonously from below. These 
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Figure 5: Local-averaged internal energy 
density at J = 0.5505. Solid lines are from 
the ordered phase, around the favored spin 
(upper one) and non- favored spin (lower one) . 
Dashed line is from the disordered phase. 
The errors are smaller than the symbol sizes. 

disordered domains increases the internal en- 
ergy density of the ordered phase and ex- 
plains why the latent heat is so small. 

Irregular boundaries of the non-favored 
clusters and disordered domains in the or- 
dered phase near the transition may be con- 
sistent with very small surface tension at the 
confined-deconfined phase boundaries found 
by a MIT-bag calculation as well as 
by lattice QCD numerical calculations |^]. 
Note, however, that the irregular boundaries 
themselves suggests inadequacy of the sharp 
spherical boundaries assumed or imposed in 
these calculations. 

The "middle-range" attractive correlation 
between the two different non-favored spins 
remains attractive down to the lowest tem- 
perature simulated, J = 0.551. Its range, 
however, decreases to mi2 = 0.07(3) at J = 
0.55075 and 0.11(3) at 0.551. This way the 
ordered phase away from the first-order tran- 
sition consolidates itself. 

It would be interesting to check if the above 
mechanism of weakening the first-order phase 
transition works in the pure-gauge QCD ther- 
modynamics. As a first step, the author gen- 
erated several hundred pure-gauge configu- 
rations on a 16^ X 2 lattice by the "APE6" 
computers in Rome, using a hybrid Monte 
Carlo algorithm written by the APE group 
1^]. From each gauge configuration a distri- 
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bution of Polyakov lines is calculated. Near 
the first-order phase transition, most of the 
Poliyakov lines are found in the neighborhood 
of the three axes with non-zero magni- 
tude, even in the disordered phase. Each 
Polyakov line is projected onto the nearest of 
the axes, resulting in distributions of the 
Potts-model spins. Correlations of the spins 
are then calculated, and found to behave in 
the same manner as their counterparts in the 
three-state Potts model. Most notably, the 
correlation between the two different non- 
favored spins in the ordered phase shows 
the same kind of "middle-range" attraction. 
This allows one to speculate that the decon- 
fined phase of QCD near the phase transition 
is susceptible to mixture of complex-shaped 
"droplets" of the confined phase . If such 
is indeed the case, perturbation calculations 
of the QCD deconfined phase would be use- 
less near the phase transition. Hence it is 
important to investigate this point in more 
detail. 

The author thanks Parallel Computing Re- 
search Facility, Fujitsu Laboratory, for com- 
puting time on their APIOOO parallel com- 
puter on which much of this work was done. 
He is also thankful for fruitful discussions 
with Norman Christ, Enzo Marinari, Giorgio 
Parisi, and Roberto Petronzio. Thanks are 
also due to the hospitality extended to him 
by the APEIOO group, especially by Federico 
Rapuano and Raffaele Tripiccione. Finally, 
he warmly thanks Andrea Donini and Ste- 
fano Antonelli for kindly allowing him to use 
their hybrid Monte Carlo simulation code for 
the APE computer. 

References 

[1] F.R. Brown et al., Phys. Rev. Lett. 
61, 2058 (1988); F.R. Brown, Nucl. 
Phys. B (Proc. Suppl.) 17, 214 (1990); 
S. Cabasino et al., Nucl. Phys. B (Proc. 
Suppl.) 17, 218 (1990); M. Fukugita, 
M. Okawa, and A. Ukawa, Nucl. Phys. 
B337, 181 (1990); Y. Iwasaki et al., 
Phys. Rev. Lett. 67, 3343 (1991). 

[2] H.W.J. Blote and R.H. Swendsen, Phys. 
Rev. Lett. 43, 799 (1979); S.J. Knak 
Jensen and O.G. Mouritsen, Phys. Rev. 



Lett. 43, 1736 (1979); M. Fukugita and 
M. Okawa, Phys. Rev. Lett. 63, 13 

(1989) ; A.J. Guttmann and I.G. Ent- 
ing, Nucl. Phys. B (Proc. Suppl.) 17, 328 

(1990) ; R.V. Gavai, ibid., 335. 

[3] M. Fukugita, M. Okawa, and A. Ukawa, 
Phys. Rev. Lett. 63, 1768 (1989). 

[4] S. Ohta, Nucl. Phys. B (Proc. Suppl.) 
26, 647 (1992). 

[5] A. Hintermann, H. Kunz, and F.Y. Wu, 
J. Stat. Phys. 19, 623 (1978). 

[6] S. Ohta, in the proceedings of the In- 
ternational Symposium on Lattice Field 
Theories "Lattice '92" in Amsterdam, 
September 1992, Nucl. Phys. B (Proc. 
Suppl.) to be published. 

[7] I. Mardor and B. Svetitsy, Phys. Rev. 
D44, 878 (1991). 

[8] J. Potvin and C. Rebbi, Nucl. Phys. B 
(Proc. Suppl.) 17, 317 (1990); K. Ru- 
mukainen, Nucl. Phys. B (Proc. Suppl.) 
26, 320 (1992); Y. Aoki and K. Kanaya, 
talk presented at the 48th annual meet- 
ing of the Physical Society of Japan, in 
Sendai, March 1993. 

[9] A. Donini and S. Antonelli, private com- 
munication. 



5 



